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abstract 

Elastic  circular  sandwich  beams  are  designed  for  minimum 
compliance  and  given  total  weight.  To  treat  the  problem  in  a  more 
realistic  maimer,  the  beams  are  regarded  as  extensible.  Examples 
are  given  for  the  optimal  design  of  circular  rings  and  semicircular 
arches  with  different  end  conditions.  The  calculated  optimal  com¬ 
pliance  is  compared  with  the  corresponding  compliance  of  a  uniform 
beam  with  identical  weight.  Finally,  the  optimal  design  with  stress 
bounds  is  also  investigated  for  the  ring  problem. 


^  The  results  presented  in  this  paper  were  obtained  in  the  course  of 
research  supported  partly  by  the  U.  S.  Army  Research  Office- 
Durham  under  Research  Grant  DA.ARO-D-3 1-61 008  and  partly  by 
the  Air  Force  Office  of  Scientific  Research,  Office  of  Aerospace 
Research,  U.  S.  Air  Force  Office  under  AFOSR  Grant  AF-AFOSR 
1226-67  with  University  of  California,  San  Diego. 


1.  INTRODUCTION 


In  the  analysis  of  the  deformation  of  the  curved  beams  such  as 
rings  and  arches,  we  usually  assume  that  the  center  line  of  the  beam 
is  inextensible  and  the  deflection  of  the  beam  is  entirely  due  to  flexure. 
Such  an  assumption  could  simplify  the  analysis  considerably  and  the 
error  introduced  by  the  assumption  is  very  small  if  the  ratio  of  the 
thickness  of  the  beam  to  the  radius  of  curvature  of  the  center  line  of 
the  beam  is  much  less  than  unity.  However,  in  the  optimal  elastic 
design  of  structures  for  given  total  weight  and  minimum  compliance, 
the  design  variable  is  proportional  to  the  strain  energy  density  [1] 
Therefore,  the  assumption  of  inextensibility  would  lead  to  a  design 
with  zero  cross  sections  at  those  locations  of  vanishing  bending  moment. 
The  fiber  stress  would  become  infinite  there  when  the  axial  force  does 
not  vanish.  To  remedy  this  deficiency  one  may  stipulate  a  minimum 
cross  section  constraint  [2]  ,  or  set  bounds  on  the  fiber  stress,  or  re¬ 
move  the  assumption  of  inextensibility. 

The  present  paper  illustrates  how  the  last  amendment  is  applied 
to  a  realistic  design.  For  simplicity,  our  investigation  is  restricted 
to  the  design  of  sandwich  beams  for  which  the  bending  stiffness  is  re¬ 
garded  as  proportional  to  the  extensional  stiffness. 
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Numbers  in  brackets  designate  References  at  end  of  paper. 


2.  BASIC  EQUATIONS 


Consider  the  in- plane  infinitesimal  deformation  of  a  circnbtr, 
curved  beam  with  symmetrical  cross  sections.  The  thickness  of  the 
beam  is  assumed  to  be  much  smaller  than  the  radius  R  of  the  center 
line  of  the  beam.  The  stress-displacement  relations  of  the  curved 
beam  can  be  derived  based  on  the  Euler- Bernoulli  assumption  [3]  . 
The  stress-displacement  relations  thus  obtained  are  too  complicated 
to  be  applied  to  practical  problems  unless  further  approximations  or 
linearizations  are  imposed.  In  the  following,  we  shall  postulate  that 
the  fiber  stress  at  any  point  in  the  beam  with  a  distance  z  from  the 
central  surface  can  be  expressed  as 


where  F  is  the  fiber  stress  resultant,  M  is  the  bending  moment, 

A  is  the  cross-sectional  area  and  I  is  the  moment  of  inertia  of  the 
cross-sectional  area  about  the  centroidal  axis.  In  fact,  Eq.  (1)  can 
be  derived  based  on  Euler- Bernoulli  assumption  if  terms  of  the  order 
of  r/R  are  neglected  as  compared  with  unity  where  r  is  the  radius 
of  gyration  of  the  cross  section  of  the  beam.  The  stress  energy  per 
unit  length  of  the  beam  can  be  expressed  by  the  following  area 
integral  over  the  cross  section: 


U 


0f»dA 


1 

2  V  AE 
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Let  the  inward  normal  distributed  load  be  p(8),  the  counter, 
clockwise  tangential  distributed  load  t(8),  the  inward  normal  displace¬ 
ment  w  and  the  counterclockwise  tangential  displacement  v.  By  con¬ 
sidering  the  equilibrium  in  moment,  radial  and  tangential  forces  acting 
on  an  infinitesimal  element  of  the  beam  as  shown  in  Fig.  1,  we  obtain 
the  following  equations  of  equilibrium: 


^  +  M'  +  p  ---  0  , 

(3) 

F'  .  +  t  =  0 

R 

(4) 

The  stress  displacement  relations  can  be  derived  from  the 
principle  of  minimum  complementary  energy  with  Eqs.  (3)  and  (4)  as 
constraint  conditions  while  the  displacements  v  and  w  are  con¬ 
sidered  to  be  prescribed.  We  obtain 


F 

AE 


(v  +  w) 


(5) 


M 

El 


(6) 


where  the  dot  represents  the  differentiation  with  respect  to  the  polar 
angle  9  .  Note  that  ( '  )  s  -  R  (  ) ' . 

Next,  let  us  consider  the  optimal  design  of  circular  curved 
beams.  The  beam  is  designed  to  minimize  total  compliance  due  to  a 
concentrated  force  of  given  magnitude  while  the  total  weight  of  the 
beam  is  prescribed.  The  necessary  and  sufficient  condition  of  local 
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optimum  is  given  in  [  1  ] .  It  is 


U/B  =  constant  (7) 

where  U  is  the  strain  (or  stress)  energy  per  unit  length  of  the  beam 
given  by  Eq.  (2)  and  B  is  the  design  variable.  In  the  case  of  optimal 
design  of  sandwich  beams  with  given  height  of  the  core,  we  shall 
choose  the  area  of  the  cover  sheets  A  as  the  design  variable.  In  this 
case,  the  stiffness  per  unit  length  of  the  beam  is  proportioxial  to  A 
and  Eq.  (7)  becomes  the  necessary  and  sufficient  condition  of  global 
optimum.  The  application  of  £q.  (7)  to  the  optinnal  design  of  straight 
sandwich  beams  is  given  in  [4]  . 

If  we  assume  that  the  curved  beam  is  inextensible,  then  the 
first  term  in  Eq.  (2)  is  dropped.  In  this  case,  the  design  variable  A 
would  vanish  whenever  the  bending  moment  M  vanishes.  Note  that 
at  any  section  with  zero  value  of  A,  the  fiber  stress  becomes  infinity. 
Thus  the  assumption  of  inextensiblity  is  no  longer  valid.  In  order  to 
treat  our  problem  in  a  more  realistic  manner,  we  shall  regard  the 
beam  as  extensible  and  use  the  expression  of  strain  energy,  Eq.  (2). 

We  shall  assume  the  cross-sectional  areas  of  the  top  and 
bottom  cover  sheets  to  be  identical.  The  core  is  regarded  as  perfectly 
soft  and  the  whole  fiber  stress  is  carried  by  the  cover  sheets.  Since 
the  thickness  of  both  cover  sheets  is  much  smaller  than  the  height  of 
the  core  T,  the  moment  of  inertia  of  the  cover  sheets  can  be  expressed 
approximately  as 
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where  P  =  T/(2R)  is  a  constant.  The  curved  beam  is  said  to  be 
slender  if  p  <<  1.  From  Eqs.  (2),  (7),  (5),  (6),  and  (8),  the 
optimality  cowlition  can  be  expressed  in  terms  of  displacement  com¬ 
ponents  as 

P*  (v  -  w +  (v  +  w)®  =  c  (9) 

where  c  is  a  positive  constant.  Note  that  Eq.  (9)  is  a  nonlinear 
differential  equation  of  v  ami  w. 

3.  OPTIMAL  DESIGN  OF  CIRCULAR  RINGS  AND  ARCHES 

A  circular  elastic  slender  sandwich  ring  of  radius  R  is  de- 
formed  by  two  equal  and  opposite  forces  P  as  shown  in  Fig.  2a.  For 
given  total  weight  of  the  ring,  we  shall  find  the  variation  of  thick¬ 
ness  of  the  cover  sheets  such  that  the  deflections  at  a  and  b  are 
minimized.  A  similar  problem  of  optimal  design  of  plastic  rings  is 
investigated  by  Prager  [2]  . 

Since  the  optimum  ring  is  symmetrical  with  respect  to  the 
diameters  0  =  0  and  0  =  ir/2,  we  need  only  consider  one-quarter  of 
the  ring  in  the  first  quadrant  0  £  0  s  ir/2.  At  0=0,  the  axial 
tension  is  P/2  and  the  transverse  shear  is  zero.  Note  that  the  bend¬ 
ing  moment  at  0  =  0  remains  unknown  and  we  shall  denote  it  by  . 
At  any  section  within  0  s  0  z  ir/2,  the  bending  moment  and  the  axial 
tension  are  respectively 


“PP 

-  —  (1  -  cos  0)  , 


F  =  —  cos  0 


After  elimination  of  M  and  F  from  Eqs.  (5),  (6)f  (10),  and 


(11),  we  obtain 


-A  ( V  +  w)  =  ®  ' 


A(w  -  V  )  =  2EP^ 


where  k  =  1  -  _ p-  is  a  constant  to  be  determined.  From  Eqs.  (9), 

PR 

(12),  and  (13),  the  design  variable  A  can  be  expressed  as 


A  PR  • 
*  '  If ' 


^  CO,®  sj  .  (14) 


Since  the  total  weight  of  the  core  is  fixed,  to  prescribe  the 
total  weight  of  the  beam  is  equivalent  to  assigning  the  total  weight  of 
the  cover  sheets.  Let  the  density  of  the  cover  sheets  be  p  .  The  total 
weight  of  the  cover  sheets  is 

W  =  4prJ  Ad0  =  c-i/®  f  (15) 


where 


-i-i/a 

f  =  J  ^  cos*  0j  d0 
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From  Eqs.  (14',  -^nd  (15)  tl>.  cros::- sectional  area  within  0  ^  9  <  ir/Z 
can  be  written  as 


A  *  (k  -  cos  9)»  +  cos»  9  (17) 

Since  the  value  of  k  is  found  to  be  nonzero,  the  value  of  A  would 
never  vanish. 

In  view  of  Eq.  (17),  the  general  solutions  of  Eqs.  (12)  and 
(13)  are  found  to  be 


w  =  c,  sin  9  +  0,.  cos  9 
1  2 


+ 


~  1  COS03  [(k-  cosp)^  +  P®  cos®  sin  (9  -  0)dp, 

ETWp  Jq 


and 


(18) 


▼  =  c  cos  9  -  c„  sin  9  +  c_ 

13  3 

.  r  /k-  COS0-  [k-  (l+p®)cosp]cos  (9-p)}  [(k  _  COB  v)® 

EWP  ’*0  ^  ■' 

+  P®  cob®^b]  . 

(19) 

The  integration  constants  c,  ,  c„,  c,  and  the  value  of  k  can  be 

a  8  3 

determined  by  the  boundary  conditions 

w  (0)  =  V  (0)  =  w  (v/2)  =  V  (v/2)  =  0  (20) 

It  is  found  that 

Cj  =  C3  =  0  ,  (21) 
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w 

I  . . .  ^ 

2  V 12,  -l/3 

C  =  r  Sind  [k>  (1+P^)C08  e]  [(k- cos8)^  +  cos^  9]  dd  . 

2  EWP  Jq 

(22) 

and  the  constant  k  satisfies  the  following  equation: 

.ir/2  .1/3 

I  (k  -  cos  d)  [(k  -  cos  d)*  +  P°  cos®  d]  dd  =  0  .  (23) 

•‘o 

Note  that  Eq.  (23)  can  also  be  obtained  directly  by  integrating  Eq.  (13) 
and  using  Eqs.  (17)  and  (20).  After  the  value  k  is  determined  from 
Eq.  (23).  f,  and  hence  A  (d),  can  be  evaluated  from  Eqs.  (16)  and 
(17)  respectively.  In  Fig.  3,  the  ratio  of  A  to  its  average  value 
is  plotted  against  d  in  solid  lines  for  p  =  0.05,  0.1,  and  0.15. 

It  is  found  that  when  P  approaches  zero  the  minimum  cross- 

I 

sectional  area  also  approaches  zero  and  its  limiting  position  is  at 
d  =  -^ .  This  is  not  surprising  because  in  this  limiting  case  the  first 
term  in  the  bracket  of  Eq.  (14)  becomes  dominant  and  the  design  is 
actually  equivalent  to  that  of  an  inextensible  ring.  One  can  easily 
show  that  for  the  optimal  design  of  inextensible  rings,  k  =  cos~^  (^/4) 
and  A(d)  is  proportional  to  |k  -  cos  d|  .  Accordingly,  we  conclude 
that  for  a  very  small  value  of  p,  the  optimal  design  of  an  inextensible 
ring  is  a  good  approximate  design  except  in  the  region  near  the  min¬ 
imum  cross  section.  We  shall  discuss  this  point  further. 

To  show  the  efficiency  of  the  optimal  design,  the  deflection 
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6  at  the  point  of  the  application  of  the  force  P  is  compared  with  that 
o 

of  a  ring  of  uniform  cross  section  and  equal  weight.  For  optimum 
rings,  we  have 

ga  6  =  -  4  f  [k- (1  +  P^)cos8]  C(k- cose)®  +  p®  cos®  0]  cos  Odp. 

2PpR  o  p 

(24) 

For  uniform  rings,  the  deflection  6^  can  be  obtained  from 
Castigliano's  theorem  as 


(25) 


These  compliances  and  their  ratio  are  plotted  against  p  in  Fig.  4. 

It  is  found  that  the  compliance  is  reduced  by  18-25%  between  p  =  0. 15 
and  0.05  due  to  the  optimtim  design. 

In  the  foregoing  discussion,  we  have  considered  the  design  of 
the  circular  rings.  The  optimal  design  of  other  types  of  the  circular 
slender  beams  can  also  be  analyzed  in  a  similar  manner.  Let  us  con¬ 
sider  the  optimal  designs  of  semicircular  arches  deformed  by  a  vertical 
force  P  at  the  crown,  as  shown  in  Figs.  2b  and  2c.  For  two  hinged 
arches,  the  bending  moment  and  the  axial  tension  can  be  expressed  in 
terms  of  the  horizontal  outward  reaction  H  at  0  =  0.  For  clamped 
arches,  the  bending  moment  and  the  axial  tension  can  be  expressed  in 
terms  of  two  redundant  quantities  and  H  at  the  end  of  the  arch. 


Tilling  appropriate  boundary  conditions  we  can  easily  derive  the  follow¬ 
ing  governing  equations  for  the  optimal  design  of  semicircular  arches: 
(a)  Two  hinged  arches 


A  =  —  (cos  e  +  h  sin  9  -  1)®  +  (cos  9  +  h  sin  0)*  J  ,  (26) 


v/2 


i/s 


g^  =  J  0  +  h  sin  0  -  1)®  +  (cos  0  +  h  sin  0)®  J  d0  ,  (27) 


w/2 

[(l+P®)(cos0  +  hsin0)  -  I][(cos0  +  hsin0  -  1)® 

0  -i/a 

+  P®  (cos  0  +  h  sin0)®]  sin9d0=O, 

(28) 

EW6  g.  ,’r/2 

pp^lf  ~  *  "p  J  [(cos  0  +  h  sin  0  -  1)  (1  -  cos  0 ) 

-  p®  cos  0  (cos0  +  h  sin0)]  [(cos0  +  hsin0  -  1)® 

-i/a 

+  P®(co8  0  +  h  sin0)®]  d0  ,  (29) 


EW6 

u 

PpR® 


2  - 


1 

w  I  +p® 


) 


+ 


w 

4 


(30) 


(b)  Clamped  arches 

A  =  ^p^  ^  0  +  h  sin  0  +  n)®  +  (cos  0  +  h  sin  0)®  J"  ,  (31) 

1  -  1/8 
gg  =  J  (cos0  +  h  sin0  +  n)®  +  (cos0+  hsin0)®J  d0  .  (32) 
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»/2 

f  C(1  ■*■?*)  (cos  0  +  h  sind)  +  n]  [(cos  0  +  h  sin0  +  n)® 

*^0  .1/3 

+  P®  (cos0  +  h  sin6)^]  sin  0d0  =  0, 

(33) 

.*/2  -i/a 

1  (cos  0  +  hsin0  +  n)  [(cos  0  +  hsin0  +  n)®+ p® (cos  0  +  hsin0)®  ]  d0  =  0, 

°  (34) 

EW6  g  ir/2 

p  -°  =  ft"  J  [(1  +  P®)  (cos  0  +  h  8in0)  +  n]  [(cos  0  +  h  sin0  +  n)® 

^  -i/a 

+  P®  (cos  0  + h  8in0)®]  cos  0d0  ,  (35) 


EW6  ,  .  -3 

ir(ir-2)[(l  +  p®)  ,r(ir+2)  -  16]/[(1  +  p®)  ir®  -  8]  . 


(36) 


2H 

In  Eqs.  (26  -  36),  h  =  -p-  and  n  =  pj^  -  1.  The  distributions  of  the 
cross-sectional  area  A  of  these  two  arches  are  shown  in  Fig.  5  and  6 
respectively.  The  comparison  of  the  compliance  of  the  optimum  arch 
with  that  of  a  corresponding  uniform  arch  is  shown  in  Figs.  7  and  8. 

By  comparison  of  our  design  with  the  design  of  an  inextensible 
beam,  it  is  found  that  the  variations  of  A  in  these  two  designs  are 
similar,  particularly  when  the  values  of  p  is  very  small.  When  the 
curved  beam  is  very  thin  the  design  is  primarily  governed  by  the  flex¬ 
ure,  except  in  the  vicinities  of  M  =  0  where  the  axial  force  dominates 
and  the  value  of  A  attains  a  local  minimum.  This  phenomenon  can  be 
visualized  by  examining  the  expressions  of  A,  Eqs.  (17),  (26),  and 
(31).  In  these  equations  the  first  term  in  the  bracket  is  due  to  flexure 
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and  the  second  is  due  to  extension.  The  first  term  is  much  larger  than 
the  second  when  p  is  small  and  6  is  not  close  to  the  value  where 
M  =  0.  When  p  increases,  the  influence  of  the  axial  force  becomes 
more  evident,  the  cross  section  becomes  more  uniform,  and  hence  the 
efficiency  of  the  optimal  design  compared  with  the  uniform  beam 
decreases. 


4.  OPTIMAL  DESIGN  OF  CIRCULAR  RINGS  WITH  AN  UPPER 
BOUND  ON  STRESSES 

Unlike  the  case  of  an  optimal  straight  beam,  the  magnitude  of 
extreme  fiber  stresses  in  an  optimal  ring  with  cross  section  as  shown 
by  Eq.  (17)  is  no  longer  uniform.  Instead,  it  can  be  written  as 


a  =  max 


^  + 

MT 

A  ^ 

21 

_  2PRpf 


-^1  +  Zp  Ik- cos0|  (k- COS0)® +COS®  pj  cos0^ 


-i/a 


In  Fig.  9,  is  plotted  against  0  for  p  =  0.  05.  It  is  found  that 

k  k 

reaches  its  maximum  value  at  0  =  cos"^  —  and  cos"^  . 

o  I  +  P  1  -  P 

and  minimum  value  at  0  =  cos”'  k  and  ir/2.  The  magnitude  of  these 
extreme  values  increases  rapidly  as  p  decreases.  Therefore,  to 
make  the  design  more  realistic  it  is  necessary  to  set  an  upper  bound 
on  the  fiber  stress. 

Therefore  the  problem  is  to  minimize  the  overall  compliance 
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C  =  a  J  eA  d0 


iinder  the  condition  of  prescribed  weight 


W  =  c  J  Ad0 


where  a  and  c  are  given  constants  and  e  is  the  strain  energy  per 

unit  volume.  In  Eqs.  (38)  and  (39)  integrals  are  taken  over  the  entire 

beam.  Since  O  is  not  greater  than  its  upper  bound  a  ,  we  have 
o  c 

a  -  ff  =  f)*  (0)  (^ 

c  o 


where  t}(0)  is  an  unknown  slack  variable  function.  By  introducing 
the  Lagrange  multipliers  and  (0)  we  have  the  following  vari. 
ational  eqviation  for  optimality, 

6  {aj  eAd0  -  xjc  jAdO  -  W)  -  Jx3(0)[a^-  T)*(e)Jde}  =  0. 


The  variations  with  respect  to  X^  ,  X^t  f).  and  A  furnish  Eqs.  (39) 


and  (40). 


as  the  Euler  equations  and 


Kv  =  0 


J  (ae  -  )  6  Ad0  +  6^  J  Xg  (a^  -  -  17*  )  d0  =  0  (43) 


In  deriving  Eq.  (43),  the  dependence  of  e  on  A  is  disregarded  because 
of  the  principle  of  minimum  potential  energy.  In  view  of  Eqs.  (40), 
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(42),  and  (43),  we  conclude  that  when  t)  =  0,  we  have 

a  =  O  (44) 

CO 

and  when  T)  0,  =  0  and  Eq.  (43)  reduces  to  Eq.  (7)  as  the 

optimality  condition.  Therefore,  for  the  optimal  design  of  rings,  if 

the  stress  bound  a  is  greater  than  the  maximum  value  of  a  ,  i.  e. , 
c  o 

then  the  stress  constraint  is  ineffective  and  the  optimality  con¬ 
ditions  Eq.  (9)  governs  the  design  of  the  entire  ring.  Hence,  (a^)  ^ 
is  an  upper  bound  of  below  which  the  stress  constraint  will  in¬ 
fluence  the  design.  On  the  other  hand,  if  is  sufficiently  small,  we 
have  either  the  case  that  the  given  weight  is  not  large  enough  for  designing 

a  ring  satisfying  the  stress  constraint  a  a  a  in  all  sections  or  the 

c  o 

case  that  the  stress  constraint  governs  the  entire  design  of  the  ring. 

In  the  latter  case,  a  =  (7  and 
o  c 

A  =  ®  J  .  (46) 

From  the  condition  of  the  given  weight,  we  find  that 

~  J  ^  0  J  d8  ■  •  (47) 

The  equation  for  k  can  be  obtained  from  Eq.  (23)  with  a  modification. 
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It  is 


r  1  n  “1 

J  (k  -  cos  0)  1 1'  -  cos  8 1  cos  ej  de 


(48) 


When  a  =  ((T  )  j  .  the  design  is  entirely  governed  by  the  stress  con- 
c  c  ^ 

straint.  When  a  <  (a  )  -  ,  the  given  weight  of  the  beam  is  not 
c  c  ^ 

sufficient  to  fulfill  the  condition  of  stress  constraint.  Finally,  when 
(T  >  ((T  )  .  .  the  stress  constraint  does  not  govern  the  design  of  the  whole 
ring  and  the  optimality  condition,  Eq.  (9),  has  to  be  taken  into  consider¬ 
ation.  Therefore,  (a  )  is  a  lower  bound  above  which  both  Eqs.  (9)  and 
(44)  have  to  be  used  to  determine  the  optimal  design.  In  Fig.  10,  the 
dimensionless  (ff  )  and  (a  ),  are  plotted  against  p.  It  is  found  that 
the  values  of  (CF^)  ^  and  (cr^)  ^  increase  rapidly  as  p  decreases. 

In  the  following,  we  are  interested  only  in  the  problem  with  o 

c 

in  the  range  (o  )  <  (j  <  (o  )  .  In  this  case,  each  of  Eqs.  (9)  and 

C  V  C  C  XI 

(44)  governs  the  design  in  some  intervals.  In  view  of  Fig.  10,  it  is 
reasonable  to  assume  that  a  =  a  in  two  intervals  8  s  8  s  9  and 
83  s  9  s  9^  where  9^  a  9^  ,  and  Eq.  (9)  holds  in  the  remaining 
intervals  of  0  s  9  s  .  Accordingly,  the  design  variable  A  will  take 
the  form  of  Eq.  (46)  in  the  former  intervals  and  the  form  of  Eq.  (14)  in 
the  remaining  intervals.  The  total  weight  of  the  cover  sheets  is  thus 

W  =  4pR  {(  J®"  +  J®*  +  )||  (k.  cos9)»  +  cos»  9]'^*  d9 

®  ®a  ®4 
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^  (K'*^  co«0]d0}  .  (49) 

®i  ®3 

Integrating  Eqs.  (12)  and  (13)  and  using  the  boundary  conditions 
Eq.  (20)  and  the  continuity  conditions 

w  (0")  =  w(0 .''’),  v(0")  =  v(0.^)  and  w(0")  =  w(0.^)  , 

i  =  1,  2,  3,  4  ,  (50) 


we  finally  obtain  the  deflection  at  0  =  —  as 
0  0  ir/2 

■'■J  (k-  COs6)  -  COB0j||-^(k-  CO8  0)® 

*  a  ^ 

+  COS*  0  I  cos  0d0 


Of  R  .  -0  .0. 


+  Q  *  +  J  *  ^  (k-  cos  0)-  cos  ®  1*^-  008  0  I 

“1-1  > 

+  cos  0  J  cos0d9^,  (51) 


®  0, 

1  3 


and 


,0,  ,0,  -■'f/2 


3  ♦ 


-  T»/» 

+  cos*  0  J  d0 


R  -0-  0 


*  “I"  (J  I  *)  "p  <*'-008  0)  P-i  |k-CO80|  +  COS0  1  ^d0  =  0, 

®1  ®3 


(52) 
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The  unknowns  c  and  k  ^an  be  solved  from  Eqs.  (49)  and  (52).  How¬ 
ever,  B  ,  6  ,  B  ,  and  d.  are  still  undetermined.  If  these  values  are 
prescribed,  the  constants  c  and  k  can  be  found  from  Eqs.  (49)  and 
(52).  Hence  the  design  is  determined  by  Eq.  (46)  in  the  intervals 
8^  s  6  s  0^  and  0^  s  0  s  0^,  and  by  Eq.  (14)  in  the  remaining  intervals 
0  £  0  s  .  The  design  of  this  type  is  optimal  among  all  designs  with 
a  =  O  in  the  intervals  0,  s  0  s£  0  and  0,  <  0  s  0  Note  that  in 

O  C  1  3  3  4 

this  design  the  value  of  A  is  not  necessarily  continuous  at  0=0^,  0^^, 
0,0.  For  the  optimal  design,  however,  the  values  of  0,  to  0,  should 
make  the  right-hand  side  of  Eq.  (51)  a  minimum  with  respect  to  all 
admissible  values.  This  criterion  finally  leads  to  a  condition  that  A 
must  be  continuous  at  0=0^,  0g,  83,  and  0^  [see  Appendix].  Hence, 
we  have 

■||c-i/='[-ig(k.cos0.)’>+  cos=0j  ‘  =  -^  [|-i  lkco8  0.|  +  cos  0^, 

i=l,2,  3,  4  .  (53) 

Equations  (49),  (52),  and  (53)  can  be  used  to  determine  c,  k,  and 
0.(i  =  1,  2,  3,  4). 

A  slightly  modified  Newton- Raphs on  iteration  method  is  used 

successfully  in  solving  these  equations.  The  broken  line  in  Fig.  3 

shows  the  value  of  A  with  the  stress  constraint  a  =  0.85  (a  )  for 

c  c  u 

P  =  0.05.  It  is  found  that  the  deviation  of  A  from  the  corresponding 


value  in  the  design  without  the  stress  constraint  is  very  small.  How¬ 
ever,  the  difference  becomes  much  larger  when  approaches  the 
value  (a  )  .  =“  0.78  (ff  )  •  The  extreme  fiber  stress  a  in  the  optimal 

beam  is  shown  by  the  broken  line  in  Fig.  9.  It  is  seen  that  the  consider¬ 
able  reduction  of  c  in  0,  s  0  s  0,  and  0-  s  0  s  0  does  not  lead 
to  a  significant  increase  of  in  the  remaining  regions. 

When  p  is  sufficiently  small,  say  less  than  0.1,  the  distance 
A0  between  two  maxima  of  in  the  design  without  stress  constraint, 
is  approximately  equal  to  2^  (Fig.  9).  In  the  design  with  stress  con¬ 
straint,  the  interval  between  two  regions  governed  by  the  stress  constraint 
is  much  smaller  than  A0  as  shown  by  the  broken  line.  Therefore,  if 
the  value  of  0  is  very  small,  it  is  reasonable  to  approximate  the  design 
by  ass\iming  that  there  is  only  one  interval,  say  0j^  s  0  s  0^,  governed 
by  the  stress  constraint. 

As  we  have  pointed  out,  the  extreme  fiber  stress  becomes  infinity 
at  the  point  of  zero  moment  in  the  optimal  inextensible  beam.  Therefore, 
to  impose  the  stress  boimd  constraint,  one  may  assume  that  there  is 
only  one  continuous  interval,  in  the  neighborhood  of  zero  moment,  gov¬ 
erned  by  the  stress  constraint.  However,  this  assumption  would  lead 
to  a  design  somewhat  equivalent  to  that  of  the  one- interval- approximation 
for  Ihe  extensible  beam  discussed  in  the  last  paragraph.  To  obtain 
more  precise  results,  one  has  to  assume  that  there  are  two  separated 
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intervals  governed  by  the  stress  constraint  even  for  the  inextensible 
beam.  This  is  due  to  the  fact  that  the  variation  of  the  extreme  fiber 


stress  9  in  the  inextensible  beam  is  similar  to  the  solid  curve  in 
o 

Fig.  9  for  any  nonzero  value  of  A  in  the  vicinity  of  zero  moment. 

5,  CONCLUDING  REMARKS 

From  the  above  investigation,  we  have  the  following  conclusions; 

(1)  A  more  realistic  optimal  design  can  be  achieved  by  con¬ 
sidering  the  extensibility  of  the  beam.  When  the  thickness- 
diameter  ratio  p  of  the  beam  is  very  small,  the  influence 
of  extensibility  to  the  design  is  restricted  to  the  vicinity  of 
minimum  moment. 

(2)  The  efficiency  of  optimal  design  defined  by  comparing  its 
compliance  with  that  of  a  uniform  beam  with  the  same  weight 
depends  on  p.  Smaller  value  of  p  leads  to  a  higher 
efficiency. 

(3)  The  stress  bound  condition  would  influence  the  design  only 
if  the  value  of  allowable  stress  is  within  a  certain  range. 

For  large  value  of  allowable  stress,  the  design  is  entirely 
governed  by  the  optimality  condition.  For  small  value  of 
allowable  stress,  the  prescribed  weight  of  the  beam  is  not 
large  enough  to  ensure  the  satisfaction  of  the  stress  bound 


condition  everywhere. 


(4)  From  the  stress  distribution  in  the  ring,  designed  without 
the  stress  bound  condition,  we  may  find  approximately  the 
regions  where  material  must  be  added  to  satisfy  the  stress 
bound  condition.  However,  in  the  regions  close  to  the 
points  of  application  of  the  load,  material  can  be  deducted 
to  fulfill  the  optimality  condition. 

(5)  By  imposing  the  condition  of  stress  bound,  the  stress  in 
the  regions  of  increased  cross  section  can  be  reduced 
considerably.  However,  the  stress  in  the  remaining 
regions  of  the  beam  only  increases  slightly. 
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APPENDIX 


A  necessary  condition  that  the  deflection  6  given  by  Eq.  (51)  has 
a  minimum  with  respect  to  fl.  (i  =  1,  2,  3,  4)  is 


/!?on  ,  ,  ^os^  ,  .  ,  ,  ,  . 

ae.  "  \  ae.  j c.k  ak  ae.  ac  ae.  “  °  ^  2. 3,  4 


(54) 


where  the  subscripts  c  and  k  indicate  the  variable  held  constant 

in  partial  differentiations,  and  ■^§-  can  be  found  from  differ- 

da.  oa . 

1  1 

entiating  Eqs.  (49)  and  (52)  with  respect  to  0. .  After  substitution. 


Eq.  (54)  yields 


e  (k  -  cos  0^)^  +  cos®  0.  J  -  [k  -  cos  0.|  +  cos  6^  Jj' 


X  |co8  (1^-  COS0.)  -  cos0jj^  e  a^®  + -^  (0^-  0^-  0^+  0^) 

+  ea3(a3  +  e  b^  )] 

+  (k- cos0^)®  +  cos®  0  J  ^  -p  |k  -  cos  0^1  +  cos  0  a^(a^ 

+  eb^)  +  ka^ag  +  e  a^  b^  j 

-  (k-cjs0.)[  eka^a^  +  e®  a^b^  +  ea^a^ 

where  the  following  notations  are  used 
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e  ■  R0  /(Eel/®)  , 


0  9  ir/2 

■  Q  J  *+  J  (k- cosfl)  (k- cos  6)*  +  cos*  0^  '  d0  , 


-i/a 


0  0  0 


/  p®i  p®s  \  r  1  “|i/® 

+  J  +J  )|_'^  (k-cos  0)*  +  cos®  0J  d0  , 


0  0,  04 


i-Va 


/  P®1  P®*  (•’’’/^  V  1  f-  1  -1"*’* 

+1  +  J  j  ^  cos®  0  (k- COS0)®  +  cos®  0  J  d0  . 


a  4 


(56) 


.0.  «0,  -w/2 


(k-  COS0)  -  cos  0  (k-  COS0)® 


a  4 


-,-1/8 

+  COS®  0  cos  0  d0  , 


.0«  .0 


L*)  -  cos  0|  +  cos  0J 


-a 


<i0  . 


®i  ®3 


bg  3  J  *  (1  -  P)  cos*  0  |k-  cos  0  I  +  cos  o”!  d( 

0 

1 

0 

+  J  •^(1+P)cos®0  |k- cos0l  +  cos  0J 


-a 


d0  , 


and  0^  <  cos"i  k  <  0^  is  assumed.  The  second  brace  of  Eq,  (55) 

can  be  easily  shown  to  be  nonvanishing  in  the  limiting 

k 

case  -»  (o^)y-  In  this  case,  0j  =  0g  =  cos’i  ^  ^ , 

03  =  04  =  cos"!  j  ,  a^  =  b^  =  bg  =  0 ,  e  =  (2)i/®  ,  and  the  second 
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brace  in  Eq.  (55)  becomes 


where  a  and  a  are  the  limiting  values  of  a  and  a.  respect- 
2  3  ^  3  3  * 

ively,  which  are  obviously  nonvanishing.  Therefore  if  we  set  the 
second  brace  in  Eq.  (55)  zero,  we  shall  not  obtain  the  optimal  design. 
Otherwise,  in  the  limiting  case,  the  expression  (57)  would  vanish. 
Hence  we  conclude  that  the  first  brace  in  Eq.  (55)  must  vanish  and 
A  is  continuous  at  8  =  6.  (i  =  1,  2,  3,  4)  for  the  optimal  design. 
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F;  Axial  tension 
M:  Bending  moment 
V;  Transverse  shear 


Forces  and  moments  on  an 
infinitesimal  arc  of  the  beam. 
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Curves  for  Optimal  Circular  Rings. 
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rches. 


^(degrtes) 

Fig.  6.  Optimal  Clamped  Semicircular  Arches. 
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Comparison  of  Optimal  Two- hinged  Arches  with  Arches  of  Uniform  Cross  Section, 


Extreme  fiber  stress  in  the  Optimal  Rings  with  and  without  Stress  Constraint. 
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